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Abstract. 

Lattice statistical mechanics, often provides a natural (holonomic) framework 
to perform singularity analysis with several complex variables that would, in 
the most general mathematical framework, be too complex, or simply could 
not be defined. In a learn-by-example approach, considering several Picard- 
Fuchs systems of two- variables "above" Calabi-Yau ODEs, associated with double 
hypergeometric series, we show that D-finite (holonomic) functions are actually 
a good framework for actually finding properly the singular manifolds. The 
singular manifolds are found to be genus-zero curves. We, then, analyse the 
singular algebraic varieties of quite important holonomic functions of lattice 
statistical mechanics, the n-fold integrals x'™'> corresponding to the n-particle 
decomposition of the magnetic susceptibility of the anisotropic square Ising model. 
In this anisotropic case, we revisit a set of so-called Nickellian singularities that 
turns out to be a two-parameter family of elliptic curves. We then find a first set 
of non-Nickellian singularities for x' 3 ' an d X > that also turns out to be rational 
or ellipic curves. We underline the fact that these singular curves depend on the 
anisotropy of the Ising model, or, equivalently, that they depend on the spectral 
parameter of model. This has important consequences on the physical nature of 
the anisotropic x'™'' s which appear to be highly composite objects. We address, 
from a birational viewpoint, the emergence of families of elliptic curves, and of 
Calabi-Yau manifolds on such problems. We also address the question of the 
singularities of non-holonomic functions with a discussion on the accumulation of 
these singular curves for the non-holonomic anisotropic full susceptibility \- 
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1. Introduction 

Singularities are known to play a crucial role in physics (particle physics p], Landau 
singularities [2j |3], critical phenomena theory, renormalization group, dynamical 
systems). They are the "backbone" of many physical phenomena, in the same way 
cohomology can be introduced in mathematics as a "skeleton" describing the most 
fundamental part of so many mathematical problem^. 

Seeking for the singular points, and/or critical manifolds of models in lattice 
statistical mechanics is a necessary preliminary step towards any serious study of 
the lattice models. If the model is Yang-Baxter integrable there is a canonical 
parametrization of the model in algebraic varieties 5 , and the critical manifolds will 
also be algebraic varieties. If one does not expect the model to be "integrable" (or 
even that the integrability of the model requires too much work to be performed), 
finding the singular manifolds of the model is an attempt to obtain, at least, one exact 
result for the model. Recalling the standard-scalar Potts model [7], it is worth 
keeping in mind that its singular manifolds (corresponding to second order phase 
transitions or first order phase transitions) are selected codimension-one algebraic 
varieties where the model is actually Yang-Baxter integrable. The crucial role played 
by the (standard-scalar) Potts model in the theory of critical phenomena, is probably 
at the origin of some "conformal theory" mainstream prejudice identifying criticality 
with integrability for two-dimensional models. 

A large number of papers [SJ HI [TU] have tried (under the assumption of a unique 
phase transition) to obtain critical, and more generally singulaiQ, manifolds of lattice 
models as algebraic varieties preserved by some (Kramers- Wannier-like) duality, thus 
providing, at least, one exact (algebraic) result for the model, and, hopefully, algebraic 
subvarieties candidates for Yang-Baxter integrability of the models. The relation 
between singular manifolds of lattice statistical models and integrability is, in fact, 
much more complex. Along this line it is worth recalling two examples. 

A first example is the sixteen vertex model which is, generically, not Yang-Baxter 
integrable, but is such that the birational symmetries of the CP\§ parameter space 
of the model are actually integrabl^, thus yielding a canonical parametrizatiorjjjl 
of the model in terms of elliptic curves This parametrization gives natural 

candidates for the singular manifolds of the model, namely the vanishing condition 
of the corresponding j-invariant (which is actually the vanishing condition of a 
homogeneous polynomial of degree 24 in the sixteen homogeneous parameters of the 
model, the polynomial being the sum a very larg jpjl number of monomials |llj). This 
codimension-one algebraic variety is, probably, not Yang-Baxter integrable. 

A second example is the triangular q-state Potts model with three-spin 
interactions on the up-pointing triangles [151 116j for which the critical manifold has 

§ And not surprisingly, cohomology is naturally introduced in the singularity theory [4]. 
X If the wording "critical" still corresponds to singular in mathematics, it tends to be associated 
with second order transitions exclusively. The singular condition for the standard-scalar q-state 
Potts model corresponds to second order phase transitions for q < 4 and first order transitions for 
q > 4. 

f We have called such models "Quasi-integrable" : they are not Yang-Baxter integrable but the 
birational symmetries of their parameter space correspond to integrable mappings 
(J A foliation of CP15 in elliptic curves. 

ft In this polynomial of degree 24 in 16 unknowns is seen as the double discriminant of a 
biquadratic. It is nothing but a hyperdeterminant [12] 13 (Schafli's hyperdeterminant 1141 of format 
2x2x2x2). It has 2894276 terms. 
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been obtained as a simple codimension-one algebraic variety [17] . This codimension- 
one algebraic variety is a remarkable selected one: it is preserved by a "huge" set 
of birational transformations [T51 Q15]. Recalling the previous "conformal theory" 
prejudice on standard-scalar g-state Potts models, it is worth mentioning that, even 
restricted to this singular codimension-one algebraic variety, the model is not Yang- 
Baxter integrable. 

People working on lattice statistical mechanics (or condensed matter theory) have 
some (lex parsimoniae) simplicity prejudice that there exists a concept of "singularities 
of a model" , the singularities of the partition function being, "of course" , the same 
as the singularities of the full susceptibility. Furthermore, they also have another 
simplicity prejudice, namely that singularity manifolds are simple sets, like points, 
(self-dual) straight lines, smooth codimension-one manifolds, the maximum complexity 
being encountered with the phase diagram of the Ashkin- Teller model [50], with 
the emergence of tricritical points [2T] 122] , forgetting the much more complex phase 
diagrams of commensurate- incommensurate models [23, 24 . This Ockham's razor's 
simplicity prejudice is clearly not shared by people working on singularity theory 
in algebraic geometry, and discrete dynamical systems [U [25l [26] (see also Arnold's 
viewpoint on singularity theory and catastrophe theory |27[). 

In fact, singular manifolds in lattice statistical mechanics (or condensed matter 
theory) have no reason to be simple codimension-one sets (or even stratified spaces). 
For lattice models of statistical mechanics, where the parameter space corresponds to 
several (complex) variables, there is a gap between a physicist's viewpoint that roughly 
amounts to seeing singular manifolds as simple mutatis-mutandis generalizations of 
singularities of one complex variable, conjecturing singular manifolds as algebraic 
varieties 8, 9 ( 10 , and the mathematician's viewpoint that is reluctant to introduce 
the concept of singular manifolds for functions of several complex variables (it is not 
clear that the functions one studies are even defined in a Zariski space). 

Singular manifolds can be well-defined in a framework that is, in fact, quite 
natural, and emerges quite often in theoretical physics, namely the holonomic 
functions [35] corresponding to n-fold integrals of a holonomic integrand (most of the 
time, in theoretical physics, the integrand is simply rational or algebraic). In Sato's T>- 
module theory [29], a holonomic system is a highly over- determined system, such that 
the solutions locally form a vector space of finite dimension (instead of the expected 
dependence on some arbitrary functions). Furthermore, holonomic functions naturally 
correspond to systems with fixed regular singularities. It is crucial to avoid movable 
singularities. For non-holonomic functions, only the ones that can be decomposed as an 
infinite sum of holonomic functions (like x, the full susceptibility of the square Ising 
model) give some hope for interesting and/or rigorous studies of their singularities. 

For one complex variable, the holonomic (or D-finite (30j [31]) functions are 
solutions of linear ODEs with polynomial coefficients in the complex variable. The 
(regular) singularities can be seen immediately as solutions of the head polynomial 
coefficient of the linear ODE, up to apparent singularities [35 . If one takes a 
representation of the linear ODE as a linear differential system, one gets rid of the 
apparent singularities, and one also sees, quite immediately, the singularities in such 
systems. More generally, for holonomic functions of several complex variables, one 
can define, and see, quite clearly, the singular manifolds of the corresponding systems 
of PDEs. In a learn-by-example approach, we will show how one can find, and see, 
these singular algebraic varieties. 



Singularities of lattice models 



4 



The paper is organized as follows. After briefly recalling the framework of the 
isotropic X s ) we will first study various examples of Picard-Fuchs systems of 
two variables associated with hypergeometric series, and generalizing some known 
Calabi-Yau ODEs 36 . We will show how the singular manifolds can be obtained 
from the holonomic systems, and from simpler asymptotic calculations. We will 
then obtain singular manifolds for quite important holonomic functions of lattice 
statistical mechanics, the n-fold integrals x^' s (corresponding to the decomposition 
of the magnetic susceptibility of the anisotropic square Ising model [37]), describing 
a set of (so-called) "Nickellian" singularities, and then getting, from a "Landau 
singularity [I] [2] approach" , a first set of other (non- Nickellian) singularities. We will 
underline the dependence of the singularity manifolds in the anisotropy of the Ising 
model. This has important consequences for understanding the mathematical, as well 
as the physical, nature of the anisotropic X ' s - The question of the accumulation of 
these singular manifolds for the anisotropic full susceptibility x, will be discussed. We 
will finally comment on the emergence of families of elliptic curves for the singularity 
manifolds, and the (birational) reason of the occurrence of Calabi-Yau manifolds on 
such problems. 

2. Holonomic functions of one complex variable: the x^'s for the 
isotropic Ising model 

Let us start with the simplest holonomic, or D- finite |30j . functions, namely the 
holonomic functions of one complex variable, by recalling important holonomic 
functions of lattice statistical mechanics, the n-fold integrals x' n ' °f the isotropic 
square lattice Ising model [3SJ 135] I3S]- These n-fold integrals correspond to the 
decomposition of the full susceptibility of the model as an infinite sum |37j of the 
n-particle contributions X ■ The singularities of these X^' s have been completely 
described and can be seen to be a very rich and complex set of points (3j 00] . In 
particular, one finds, in some well-suited variable k, which is the modulus of the 
elliptic function parametrizing the two-dimensional Ising model, that the unit circle 
\k\ = 1 will be a natural boundary for the full susceptibility x °f the Ising model |40) . 
The singularities of the x^'s accumulate on the unit circle. This is the reason why 
we have this unit circle natural boundary jSHJ HOI S3 EH S3] for the full magnetic 
susceptibility x- Singularities also accumulate inside the unit circle (see Figures 1, 2, 

3, 4 of ED]), probably becoming an infinite set of points dense in the open disk \k\ < 1. 
They also accumulate outside the unit /c-circle |fc| > 1, probably becoming another 
infinite set of points also dense outside the unit circle \k\ > 1. This accumulation of 
singular points of the linear ODEs of the x^ ' s is thus (probably) dense in the whole 
k-complex plane. In other words, we do have an infinite set of singularities dense in 
the whole k-complex plane. This seems to confirm the mathematician's reluctance 
to consider singular manifolds of functions of several complex variables that are not 
holonomic: even in the very simple case of one complex variable, we already seem 
to encounter serious troubles. The full susceptibility x, which is an infinite sum [37] 
of these X s ; does not even seem to be defined in a Zariski space. Recalling these 
results [40] , the common wisdom identifying the singularities of the partition function 
and the singularities of the full susceptibility is no longer obvious. 

There is, however, an important subtlety here: these singularities are singularities 
of the linear ODEs of the X s > but not of the (series expansions of the) x^'s given 
by holonomic n-fold integrals. When one considers the fc-series expansions for the 
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's, one finds out that the singularities inside the unit circle in the open disk 
< 1, are not singularities of these series [40] . This is a quite non-trivial result. 
This is also the case for the fc-series expansion for the full susceptibility \ which is 
the infinite sum of the X s - F° r the full susceptibility x, the accumulation of X s 
singularities on the unit circle makes this unit circle a natural boundary [40] . Switching 
from high-temperature series expansions to low-temperature series, we have a similar 
result for \k\ > 1. We thus have a quite drastic difference between the singularities 
of the n-fold integrals x^ ? which are solutions of linear ODEs (they are D- finite or 
holonomic, sec below), and the full susceptibility x which is not solution of a linear 
ODE (it is not holonomic). 

Before generalizing to several complex variables with the case of the x (n) 'sforthe 
anisotropic square Ising model with two complex variables, let us consider, in a learn- 
by-example approach, simple Picard-Fuchs systems associated with hypergeometric 
series of two complex variables. 



3. A first simple Picard-Fuchs system with two variables 

Let us consider the double hypergeometric series, symmetric in x and y 

n— rn—0 

= £ ^Jr ■3F 2 ([n + l, n+i n+~], [1, l];27y)-x n (2) 

n=0 

= 1 +6-(x + y) + (90-{x 2 +y 2 ) + 720-xy) 

+ (1680 • (x 3 + y 3 ) + 45360 • x y ■ (x + y)) + (34650 • (a; 4 + y 4 ) 
+ 2217600 • x y ■ (x 2 + y 2 ) + 7484400 • x 2 y 2 ) + ■■■ 
This series reduces, when y — x, to 



\- r (3w)! ^A* 

L(n!) 3 ^ \k 
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1 + 12 ■ x + 900 • x 2 + 94080 ■ x 3 



+ 11988900 -x 4 + 1704214512 ■ x 5 + 260453217024 ■ x e +■••, (3) 

which is the solution analytic at x = of the order-four Calabi-Yau operator f2 
introduced by Batyrev and van Straten (section 7.1 of [3S], see also the ODE number 
15 in 01]) 

fl = 9 4 - 3x- (7 6 2 + 7 9 + 2) • (3 6* + 1) • (3 + 2) 

- 72a; 2 ■ (3 6 + 5) • (3 6 + 4) • (3 6 + 2) • (3 6 + 1), (4) 

where: 6 = x ■ ——. 

dx 

The double hypergeometric series ([1]) is the unique analytical (in x and y) 
solution of the Picard-Fuchs system corresponding to the two partial linear differential 
operators: 

Cl x = 9 3 x - x ■ (3 6 X + 3 6 y + 1) ■ (3 6 X + 3 6 y + 2) • (3 6 X + 3 6 y + 3), 

n y = 9 3 - y(36 x +38 y + 1)-{36 X + 39 y + 2)-(36 x +36 y + 3), (5) 



Singularities of lattice models 6 

where: 9 X = x ■ ^ 6y = V m J^- 

The other formal series solutions of ([5]), around (x, y) = (0,0), have the form 

H (x,y)-HxT -Hy) m + ••• (6) 

where the maximum value reached by n and m is 2. They read for instance: 

H (x, y) ■ ln(ir) + H x (x, y), H (x, y) ■ ln(y) + Hx(y, x), 

H (x, y) ■ ln(x) ■ ln(y) + H x (y, x) ■ ln(x) + H^x, y) ■ \n(y) + H 3 (x, y), 

It is crucial to note that the dimension of the space spanned by these formal series 
is finite. In the case of the Picard-Fuchs system ([5]), the number of solutions (i.e. 
dimension) is nine. These nine formal solutions are given in | Appendix A| The double 
series analytic in x and y, Hj(x, y) are either symmetric like Hq(x, y), H^(x, y), or 
are not symmetric like Hi(x, y). 

Such holonomic systems are also called D- finite [301 I3T] . for that reason: 
remarkably, they have a finite number of independent solutions, in contrast with 
generic systems of PDEs that have, generically, an infinite number of solutions. 
Systems of PDEs can also have no solution at all. Generically the compatibility 
of the two operators Q x and VL y , requires some (slightly tedious) differential algebra 
calculations. 

One can also see the system (J5J) as a (two-dimensional) recursion: 

(n + l) 3 • Cn+l,m = b(fl, m) ■ C n ,m, 

(m + l) 3 • c„ !m+ i = b(n, m) ■ c„, m , where: 

b(n,m) = (3(ra + m) + 1) -(3(71 + 771) + 2)- (3 (n + m) + 3), (7) 

Here, the compatibility between the two partial differential operators Q x and tt y is 
easier to see at this (double) recursion level. Introducing 

b(n,m) c n+ i m b(n, m) c„, m+ i 
ai(n, m) = —3 = , a 2 {n, m) = — T = — , 

we have the identity: 

0:2(71, m) ■ Oi(n, 771+1) = ai(n, m) ■ a 2 (n + 1, m), (8) 

which, from a recursion viewpoint, actually corresponds to the compatibility between 
the two partial linear differential operators f2 x and fl y . 

The discriminant of the two-parameter family of Calabi-Yau 3-folds read^l (see 
Prop. 7.2.1 of [5B]): 

(x + y) 3 - 3- (x 2 -7xy + y 2 ) + 3 • (x + y) -1, (9) 

or, (without performing the (x, y) —> (x/27, y/27) rescaling mentioned in |36j ) : 

A= 19683(x + y) 3 - 2187 • (y 2 + x 2 - 7 x y) +81- (x + y) -1. (10) 

This expression can easily be obtained as the resultant [12] in A (or equivalently in 
B) of the two (very simple) homogeneous binary cubics |36j : 

27 x ■ (A + B) 3 ~ A 3 = 0, 27y- (A + B) 3 - B 3 = 0. (11) 

I Note a misprint in Prop. 7.2.1 of l.'iiil : (x + y) must be changed into 3 ■ (x + y) . 
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3.1. Singular manifolds 

What are the singularities of the double hypergeometric series like (P), and how do 
they compare with the singularities of the Picard-Fuchs system ([5]) , assuming that the 
notion of singularities of such PDEs systems is well-defined ? 

From a mathematical viewpoint, when introducing some "canonical" system, 
equivalent to the Picard-Fuchs system, one should "in principle" be able to see 
the singularities as simple poles of this equivalent system. Unfortunately, to our 
knowledge, the implementation of such procedure is available as formal calculation 
tools is still in development [32] (see also [33l [34] V 

A physicist's down-to-earth approach amounts to reducing the double 
hypergeometric series, like ([1}, to scries in one (complex) variable imposing some 
relation between x and y, compatible with the (x, y) = (0, 0) origin of the double 
series. Imposing, for example, y — cx (c = 2, 3, ...), or y — ci 2 , one gets a series 
in one (complex) variable x and, then, in the second step, finds the corresponding 
linear ODE annihilating this series. The head polynomial of the corresponding 
linear differential operator gives (after getting rid of the apparent singularities) the 
singularities of these linear differential operators. An "accumulation" of such results 
enables to see that the singularities are always on the (genus-zero) algebraic curve 
S(x, y) = 0, where 

S(x,y) = 3 9 -(x + y) 3 - 3 7 • (y 2 + x 2 - 7x y) + 3 4 ■ (x + y) -1, (12) 
which is nothing but the discriminant (|10[) of the two-parameters family of Calabi- 
Yau 3-folds previously mentioned |36) . Remarkably, but not surprisingly, the singular 
variety has an interpretation as a fundamental projective invariant |12] . 

The (genus-zero) singular curve (|12j) can be parametrized by 
/l \ 3 /l \ 3 

x = (e + u ) ' v = (e - u ) • (13) 

or 

*<«>= (^>) 3 < »<«) - (e^r -*(;)■ 

where the Atkin-Lehner-like involution u <H> 1/u could suggest a modular curve 
interpretation of (TT2"j) . 

The accumulation of calculations is quite tedious compared to the simplicity of 
the final result (fi"2l) . It is far from obvious that ([121) is the singularity manifold of the 
double series ([1]), or the singularity manifold of the Picard-Fuchs system ([5]). Let us 
find a Picard-Fuchs system for which it will become crystal clear that (|12p is actually 
the singularity manifold of the system. 



3.2. Other representations as PDE systems 

In fact, the Picard-Fuchs partial differential system ([5]) can be recast into a system 
of two differential equations, each one being a linear ODE on only one variable. 
We considci[3 a linear combination of Q x , tt y and their derivatives, and cancel the 
coefficients in front of the undesired derivatives. We obtain the following form 

9 9 

fix = £ P n( x > y) ■ D *> & y = Q-(x, y) ■ D n y , (15) 

n=0 n=0 

f For our purpose, we did not use the Groebner basis approach (use the pdsolve command on the 
system of equations obtained from the Rosenfeld-Groebner command in Maple). 
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d d 
where: D x = -, D y = -, 

where P n (x, y) and Q n (x, y) are polynomials of the two variables x and y. The 
partial differential operator Cl x can be seen as a linear differential operator in x 
depending on a parameter y (and similarly f2 y as a linear differential operator in y 
depending on a parameter x). The polynomials P n (x, y) appearing in Q x will not 
be given here. For Pg(x, y) the monomial of highest degree in x and y is a; 15 y 9 (see 
(fTBj) and (|B.2|) in |Appcndix B[ ), and, for P%(x, y), ■■■ , Po(x, y), it reads, respectively, 
x 14 y 9 , x 13 y 9 , x 12 y 9 , x n y 9 , x 10 y 9 , x 9 y 9 , x 8 y s , x 7 y 7 , x 6 y 6 . 

There is a "price to pay" to recast the Picard-Fuchs partial linear differential 
system (JSJ) into a system like (TTSj) . The partial linear differential operators Cl x and 
Q y are much more involved than operators £l x and Q y in ([5]), and of higher order in 
D x or D y . The operator Q x (resp. Q y ) is of order nine with respect to D x (resp. 
D y ), in agreement with the previously mentioned finite set (IA.2I) of nine formal series 
solutions of the Picard-Fuchs D-finite system ([5]). We have checked that these nine 
formal solutions (|A.2[) are indeed solutions of Q x (resp. Cl y ). 

As a consequence of the exact symmetry interchange x <-> y of (Q]) , the partial 
differential operator Q y is nothing but operator £l x , where x and y are permuted. 
Not surprisingly, the head polynomials in (1151) have the form 

Pd(x, y) = x 6 ■ Vd(x, y) ■ S(x, y), Qg(x, y) = y 6 ■ P 9 (y, x) ■ S(x, y), (16) 

where Vg(x, y) is a polynomial of x and y, corresponding to the apparent singularities 
of the ( y-dependent) linear differential operator Q x . The expression of Vg(x, y) is 
given in | Appendix B| 

3.3. Operator factorizations 

One can actually go further in the analysis of these order-nine operators. The order- 
nine partial linear differential operator Q Xl in fact, factorizes in three order-one 
operators, and an order-six operator: 

q = (d - dln(ri(^ V)) ^ ( D d\n(f 2 (x, y)) 



dx ) V dx 
dln(r 3 (x, y)) 



(p x 



) • L 6 (x, y), (17) 



dx 

where the order-six operator L$(x, y) reads 

i 6 

L 6 (x, y) = — • Y, Pn[x, y) ■ D n x , (18) 

Pe(x, y) n=Q 

and where ri(x, y), f 2 {x, y) and f 3 (x, y) are rationals functions of x and y, while 
Pe{x, y) has simple factorizations: 



r i(x, y) = -g — —, r , r 2 {x, y) - 



x & -S(x,y)-qi " ' " a; 5 • S{x, y) ■ q 2 ' 

h(x, y) = — — j7 g , 2 , r, Pe{x, y) = x 4 ■ S(x, y) ■ V 6 (x, y), (19) 

x 4 • S[x, y) ■ Pe[x, y) 

where Vq(x, y), Vq(x, y), qi, q 2 , are polynomials of x and y given in | Appendix B| 
Not surprisingly the {x, y)-asymmetric polynomials Ve{x, y) and Vg(x, y) correspond 
respectively to apparent singularities of the order-six and order-nine operators 
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L 6 (x, y) and Q x . The polynomials p n {x, y) appearing in L e (x, y) will not be given 
here. For Pq(x, y) the monomial of highest degree in x and y is x 13 y 9 (see (|19p 
and (|B.3[) in | Appendix B[ ), and, for Ps(x, y), ■■■ ,Po(x, y), it reads, respectively, 



x 13 y 9 , x 12 y 9 , 



x 11 y 9 , 



x 10 y 9 , x 9 y 9 , 



X"y", x 7 y 7 . 



Do note that the critical exponents of this order-six operator Lq(x, y) are 
independent of y. For instance at x = the indicial polynomial reads P(r) = 
r 3 • (r — l) 3 . More remarkably, on the singular variety S(x, y) = 0, the critical 
exponents of Lq(x, y) are also independent of y. The indicial polynomial, at 
S(x, y) — 0, reads P(r) = r ■ (r — l) 2 • (r — 2) • (r — 3) • (r — 4). The singular 
behaviour at <S(a:, y) = is thus logarithmic. The wronskians of this order-six linear 
differential operator Lq(x, y), and of the order-nine operator Q x are rational functions 
of x and y, which read respectively: 

Pg(x, y) \ V 9 (x, y) 



W 



\Lq{x, y) 



W 



(20) 



x 12 ■ S(x, y) 4 ' " V"" x 7 x 27 ■ S(x, y) 7 ' 

In fact, the operator Lq(x, y) is not only Fuchsian with rational exponents and 
rational wronskian, it is actually globally nilpotent for any rational values of y. The 
p-curvature of this globally nilpotent order-six operator, is a nilpotent 6x6 matrix 
which can be put into the following Jordan fornip, not only for any rational value of 
y, but, actually, for any y being an algebraic number. 



C = 



where: 



C 4 = 0. 



(21) 



1 

1 
1 
1 


Furthermore Lq(x, y) is such that its exterior square is of order fourteen (D x 4 ) 
instead of the order fifteen one should expect generically for an order-six irreducible 
operator. This remarkable property is related to the fact that Lq(x, y) is homomorphic 
to its (formal) adjoint, with an order-two intertwinner differential operator l2(x, y) 

Lq(x, y) ■ h(x, y) — adjoint(Ji(x, y)) ■ adjoint(Le(x, y)), where: 



h{x,y) = 3 6 



27x + 27y +2 



D 2 X + Ri(x,y) ■ D x + Ro(x,y), (22) 



S(x, y) 

where R\{x, y) and R,2{x, y) are rational functions of x and y. 

One can check that the double (a;, y)-symmetric series ([1]), solution of the order- 
nine operator f2 x , is, in fact, annihilated by the order-six linear differential operator 
Lq(x, y) and, thus (by x y symmetry) by the other order-six operator 



L(>(y, x) 



Pe(y, x) 



6 

E 

n=0 



Pn(y, x) 



(23) 



At this step, we should recall that our purpose is to get the singularities of the 
system ([5]) and not to obtain an equivalent system for ([5]). Generically, systems of 



f Of characteristic polynomial P(A) = A 6 and of minimal polynomial P m (A) 
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linear PDEs cannot be strictly recast into a form like (|15[) . even for D-finite systeniijf| 
The two order-six operators L 6 (x, y) and L e (jj, x) form a PDE system that is not 
equivalent (in the sense of equivalence of systems) to the Picard-Fuchs system (0). 
However, and as far as the double series Ha(x,y) is concerned, the three systems 
(Cl x , fly), (fl x , fly), or (Lq(x, y), L§(y, x)), can alternatively be considered. 

Remark: Recovering the Calabi-Yau order-four ODE ((4|) from the y = x 
limit of the Picard-Fuchs system ([5]), or (|15j) . is not straightforward (as one could 
naively imagine). Within the (down-to-earth) approach which amounts, for instance, 
to restricting to the straight lines y = c ■ x, where c is a constant, and finding the 
linear differential operator in x, one obtains an order-six linear differential operator 
with coefficients that are polynomials in x, as well as in the constant c. One can, 
then, take the c — > 1 limit and actually recover the Calabi-Yau order- four ODE (j4]). 
These calculations are displayed in |Appendix C| The (genus-zero) singular curve (TT2"]) 

(1 - 108- (x + y)) • (2 + 27- (x + y)) 2 + 3 9 • (x - yf = 0, (24) 

reduces, in the y — x limit, to (1 — 216 a;) • (1 — 27 x) 2 — 0, namely the singularities 
corresponding to the order- four Calabi-Yau ODE (UJ). 

4. More Picard-Fuchs systems with two variables 

Similar calculations can be performed with double hypergeometric series generalizing 
the analytic solution of another Calabi-Yau order- four ODE (see |Appendix D| below). 
One can perform exactly the same calculations mutatis mutandis. 

4-1- More Picard-Fuchs system with two variables 

Let us, first, consider a two-variables Picard-Fuchs system "above" another Calabi- 
Yau ODE [35] (see the ODE number 16 in appendix A of [4"5]). corresponding to the 
following (x, y)-symmetric series with binomial coefficients: 

y, ~ /2n + 2m\ fn + m\ 2 /2n\ /2m\ „ m 
^To ™~o \ n + m )\ n ) \n ) \m J V 

m = v 7 

= 1 + 4 (x + y) + (36 (a; 2 + y 2 ) + 96 x y) + [2160 {x 2 y + x y 2 ) + 400 {x 3 + y 3 )} 
+ [4900 (x 4 + y 4 ) + 44800 {x y 3 + x 3 y) + 90720 x 2 y 2 ] + • • • (25) 
This hypergeometric double series is solution of the Picard-Fuchs system of PDEs 

n x = e 3 x -Ax- (20 x + i) (e x + e y + i) (2^ + 2^ + 1), 

fly = 6 3 y -Ay (2fl w + l) (d x + e y + l)(26 x + 26 y + l). (26) 

| Non-holonomic systems cannot be recast into a form like J 1 5 1) - This is the case, for instance, of 
the system of linear operators (fl x , Qy) = (D^., D x D y ), which has an infinite number of solutions, 
namely c ■ x + f(y) where f(y) is an arbitrary function of y. 

§ For instance, the solutions of the D-finite system (£l x ,Q y ) = (D^ —yDy, D x D y ) are solutions of 
the D-finite system (il x , Cl y ) = {D x , yDy + Dy), but this last D-finite system has more solutions. 
One needs additional operators, to have a system equivalence. 
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In the y = x limit, this series reduces to the series 




= 1 +8x + 168a; 2 + 5120 x 3 + 190120 x 4 + 7939008 a; 5 (27) 
+ 357713664 a; 6 + 16993726464 x 7 + 839358285480 x 8 + ■■■ 

annihilated by the order-four Calabi-Yau operator: 

9 4 -Ax- (5(9 2 + 56» +2) • (29 + l) 2 

+ 64a; 2 • (29 + 3) • (26> + 1) • (29 + 2) 2 . (28) 

The recast of the PDE system for the double series (f25j) into the form (fl5|) . gives 
two (x, y)-symmetric linear differential operators of order nine. The singularities of 
the two order-nine linear differential operators are respectively a; • (1 — 16 a;) = 
and y • (1 — 16 y) = Q together with the quadratic condition: 

S 2 (x,y) = 2 s -(x~y) 2 - 2 5 • (y + x) +1 = 0, (29) 

which has the simple rational parametrization 

The singularities S 2 (x, y) = are, here also, logarithmic, the local exponents being 
0, 1, 1, 2, 3, ••• , 7. 

These two (x, y)-symmetric order-nine operators also factorize in exactly, the 
same way as (fl"7]) . in three order-one operators and an order-six operator like (fT8|). 
The exterior square of this order-six operator is also of order fourteen (instead of the 
order fifteen one expects for a generic irreducible order-six operator), and, again, this 
order-six operator is homomorphic to its adjoint with a relation similar to (|2"2")l . the 
head coefficient in the order- two intertwinner being replaced by 2 8 (16a; — 16 y + 
3)/52(a;, y)/(16x — l)/a; 2 . We also have relations similar to (f20| for the various 
wronskians. 



4-2. Another Picard-Fuchs system above the Calabi-Yau operator 

Note that the Picard-Fuchs system of two variables "above" the Calabi-Yau operator 
(|28| is not unique. Other (x, y)-symmetric series reduce to the series (|27j) annihilated 
by (|28|) . for instance, the double series expansion: 



oo oo 



V- V- ,-.,<■ m. (1/2)^- (l/2^- (1/ 2),,,. n 
n = 0m = o (lf n+m -n!m 



^ ^ 64 n+m • v ' ,n n ] l ' ' • x n y m (30) 

(4)r- 



^ — ' v m 

m = 



Q£) - (31 » 



:( ^ \, ~, ~ + m], [m + l,m + l,m + 1]; 64a;) ■ f(>4//)' 
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= 1 + 4- (y + x) +3- [27- (a; 2 + y 2 ) + 2 ■ x ■ y] 
+ 20 • (y + x) ■ [125 • {x 2 + y 2 ) - 122 • x ■ y] 

+ 35/16 • [42875 • (x 4 + y 4 ) + 162 • x 2 ■ y 2 + 500 • x y ■ (x 2 + y 2 )] 
+ 63/4 • (y + x) ■ [250047 • (x 4 + y 4 ) - 248332 • x y ■ (x 2 + y 2 ) 
+ 248602 ■ x 2 y 2 ] -\ 

where (a) n is the usual Pochhammer symbol. This series can be found in Guttmann 
and Glasser [53] as a lattice Green function. It can also be seen as the expansion of a 
Kampe de Feriet function [46l H3 |48j [49] (see |Appendix D[ ): 

F S;o$ ([5]. 5. 5]. [5. ^ 5]; [1,1,1],-,-; 64a;, 64 y ) . (32) 

The double series ([30]) is not a series with integer coefficients but it can be 
recast into a series with integer coefficients if one performs the simple rescaling 
(x, y) — > (4s, Ay). One obtains: 

1 + 16 • (x + y) + [1296 • (x 2 +y 2 ) + 96-x-y] 

+ 1280 -(y + x)- [125 • (x 2 + y 2 ) - 122 • x ■ y] (33) 
+ [24010000 • (x 4 + y 4 ) + 280000 • x y ■ (x 2 + y 2 ) + 90720 • x 2 ■ y 2 ] 
+ 16128 -(y + x)- [250047 • (x 4 + y 4 ) - 248332 • x y ■ (x 2 + y 2 ) 
+ 248602 • x 2 ■ y 2 ] + ■■■ 

The recast of the PDE system for the double series into the form (fT5|) gives 
two (cc, y)-symmetric linear differential operators, now, of order thirteen. 

The singular varieties of the two order thirteen operators Q x and Q y are 
respectiveljol x ■ (x — y) ■ (1 — 64.x) = and y ■ (x — y) ■ (1 — 64 y) = 0, 
together with a (x, y)-symmetric genus-zero biquadratic which reads: 

S 2 (x, y) = 2 12 - x 2 y 2 - 2 7 ■ xy (y + x) + (x ~ y) 2 = 0. (34) 

The local exponents at the singularities of the order thirteen partial linear differential 
operators are independent of y (respectively x). 

This genus zero curve (1341) has the rational parametrization (well-suited for series 
expansions near (x, y) — (0, 0)) 

x(t) = u 2 , y(t) = (-^-Y, (35) 

M + 8u/ 

or the rational parametrization 

the Atkin-Lehner-like involution u <-> 1/u suggesting a modular curve 
interpretation of (f3~4"|) . 

Note that the two singular varieties 5*2 (x, y) and ^(x, 2/) (see (1291) ). are related 
by a simple involution: 

S 2 (x,y) = 2". x 2 y 2 - 5 2 (^, ^-). (37) 

J Such series are called globally bounded I50| . 

Note that the limit y = x of the Picard-Fuchs systems associated with (130 D . is actually a singular 
limit. 
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We thus see that the various Picard-Fuchs systems "above" a given Calabi-Yau 
ODE, (i.e. reducing, when one takes the "diagonal" y = x, to the same Calabi- 
Yau ODE), do not have necessarily the same singular manifolds, even if these various 
singular manifolds must reduce to the same singular points in the y — x limit. Since 
the singular variety ([34]) contains the origin (x, y) = (0, 0), it is easy to find, using 
the parametrization (|35|) . a linear differential ODE satisfied by (|30|) when restrictedj] 
to the singular variety (|34|) (see (|D.13j) in |Appendix D[ ). This cannot be done for ([29]) 
which does not contain the origin {x, y) — (0, 0). 

Breaking the [x, y)-symmetry in (|30p , by resumming the series as (|31[) . 
corresponds to the viewpoint of seeing Kampe-de-Feriet functions of several complex 
variables as straight generalization of hypergeometric functions [46] 1471 148] 149] . The 
^-singularities in each of the (transcendental) 4F3 coefficients of the y-expansion ([31]) 
are only the well-known af = 0,a; = l,af = oo singularities of hypergeometric 
functions (here x = 1 becomes x = 1/64), and are, of course, drastically different 
from the singular variety (|34|) for the double series ([30]) . 

The results for (f30|) . can be generalized to more general (Kampe de Feriet) double 
series depending on several parameters. 



K(x,y) = EE (a y/J'^ ' / i/ " i:]S) 

n = Q m = Wm+n n - " l - 



where (ct) n is the usual Pochhammer symbol. The same calculations as before show 
that their singular curves do not depend on the parameters. These calculations for 
([38]) are displayed in | Appendix D[ 

4-3. Picard-Fuchs systems with more than two variables "above" the Calabi-Yau 
operator (28\) . 



For heuristic reasons, we restricted to two variables but one can find many Picard- 
Fuchs systems, with more than two complex variables, "above" a given Calabi-Yau 
ODE like ([28]) . For instance, the series ([27]) of the Calabi-Yau operator ([28]) can also 
be written as the x = y = z = t subcase of the (hypergeometric) series of four 
complex variables |36j : 

'2ti + k + i + m)\ ( (3 + k + i + m y. yi (39) 

, j + k + l + m J V j\k\l\m\ ) \ X V • 

The general term being hypergeometric, one obtains directly a system of four PDEs, 
from which we build a linear ODE in the variable x, with y, z and t as "parameters". 
Once one has series with four variables, and systems of PDEs with four variables, one 
can take many limits in order to reduce to two variables. 

For instance, if one restricts the previous series to y = z = t, one gets a series 
of two variables (which will of course reduce, for y — x, to the series ([2"T]) of the 
Calabi-Yau operator (|28[)). but is no longer symmetric in x and y. The series can be 
written as: 

(2N\ 

N=0 



EL)' 3 F 2 ([-N,-N,1/2},[1,1};4). 2 F 1 ([N + 1, N + 1/2], [l];4x)-y N 



f See also the notion of Fuchsian system of linear partial differential equations along a submanif, 
(see [51) . in particular paragraph 6). 

§ The parameters of the hypergeometric functions become linear differential operators 48 49 . 
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= E(/v) ■ ^F 2 ([-N,-N,-N}, [1,1/2 -N]; 1/4) 
x 2^(^+1, TV + 1/2], [1]; 4x)-y N 

= E E u j?' • 3F 2 ([-m,-™,l/2], [1,1], 4) (40) 

n=0 m=0 ^ ' 

The corresponding system of PDEs reads: 
n x = e 2 x - 2x- (6 x + 9 y + l)(20 x + 29y + l), 

n y = 9 y - 2y {109 2 y + 106 y + 3){9 x +9 y + l)(29 x + 29 y + l) (41) 

+ 36 y 2 ■ {26 x + 29 y + 3) (29 x + 2e y + 1) {6 X + e y + 2){8 X + 9 y + 1). 

Again, one can recast this system into a form like (H"5j) . i.e. two linear differential 
operators £l x and Q y in the variable x (resp. y), both of order eight, each one with 
the same singular variety which is the union of the two genus-zero algebraic curves: 

16a; 2 - 8- {Ay + 1) • x + {Ay -if = 0, and: 

16 a; 2 - 8- (36 y + 1) • x + (36 y - l) 2 = 0. (42) 

These two order-eight operators both factorise in a similar way as p7|) but, this time, 
in the product of two order-one and one order-six operator. These two order-six 
operators rightdividing respectively fl x and Cl y are not related by a (a;, j/)-symmetry, 
because the Picard-Fuchs system (1411) is not (x, y)-symmetric. Again these two order- 
six operators are such that their exterior square are of order fourteen instead of the 
order fifteen one can expect for the exterior square of a generic irreducible order-six 
operator. Furthermore one has, again, that these order-six operators are homomorphic 
to their adjoint, the intertwinner being of order two (see (f22l ). 

More examples of Picard-Fuchs system with two variables "above" Calabi-Yau 
ODEs are sketched in |Appcndix E| their corresponding (simple) singular varieties 
being also given. 



5. Singular manifolds for hypergeometric series of several complex 
variables 

All these singular varieties (TT2l . (|24l) . (|34]l (as well as similar ones, (|E.3[) . (|E.7I) . given 
in | Appendix E I can, in fact, be easily obtained from very simple calculations when one 



remarks that the previous double series are hypergeometric series of several complex 
variables. The calculations, corresponding to the Horn's convergence theorem, are 
similar to the ones for Horn functions and Horn systems [551 [531 (Ml [55J. A very 
important property is the fact that the region of convergence for hypergeometric series 
does not depend on the parameters [56] . 

Let us denote the coefficients of ([T|), by c n ,m 

(3m + 3n)! 



the successive ratio of c n , m in the two "directions" reads respectively 
c n ,m (n+1) 3 c„, m (m + 1) 3 



(43) 



Cn+1,TO b(n, m) c n , m+ i b{n, m) 
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where the product 6(n, to) is given by ([7]). In the n and to large limits these two 
ratios behave respectively like 

77 777 

-XYro, m) = ; 7-=-, and 5^fn, ml = ; rr, (45) 

V ' 7 27 (m + n) 3 ' V ' ' 27 (m + 77) 3 ' 1 ' 

where one remarks that X(n, to) and Y(n, to) depend only of the ratio 72/777. The 
curve rationally parametrized by (a;, y) = (X(n, to), K(n, m)) can easily be obtained 
performing a resultant (elimination of to or 77 or the ratio 72/777) and one recovers, in 
a very simple way the singular manifold (fT2"|) . One notes that ([4"5"f is nothing but the 
previous binary cubics ([TTj) yielding (fTTj|) . the discriminant of a two-parameters family 
of Calabi-Yau 3-folds. 

We can perform similar calculations for the hypergeometric series (|25l) . the ratio 
of the c„ im 's also read (|44|) . the product 6(71, to) being now given by 

o(?2, TO) = 2- (2 72 + 2m + 1) (2 72 + 2 772 + 2) (2 72 + 1). (46) 

In the 77 and to large limit, this gives the rational parametrization of the singular 
variety (1291) . namely (x, y) = (X(n, m), Y(n, to)), with: 

2 2 
77 TO 

m ) = TJ7 1 — \2' and m ) = ift ^ , — ^2- ( 47 ) 

lO (777 + 72) z lO (777 + 72) z 

For the hypergeometric series (|30)) . the ratio of the c„ iTO 's read respectively 

(77 + TO + 1) 3 (72 + 1) (72 + TO + l) 3 (TO + 1) 

4- (2 72 + 1) 3 (2 72 + 2 TO +1)' 4- (2 777 + 1)3 (2 77 + 2 TO +1)' 

In the 72 and to large limits this gives the rational parametrization of the singular 
variety ([34]), namely (x, y) — (X(n, m), Y{n, to)), with: 

X(n,m) = ^±#, and F(n, m) = i^±#. (48) 

d4 77^ 04 to z 

Finally for other hypergeometric series (IE.2[) . (|E.6[) . given in |Appendix E| 
similar calculations also give rational parametrizations of the corresponding genus- 
zero singular curves (|E.3[) and (|E.7I) . 

For instance the successive ratio of c„. m 's for (|E.6|) read respectively 

h T' T7 V' where: ( 4Q ) 

0(72, TO) 0(72, 772) 

b(n, m) = (2 77 + TO + 1) (2 72 + TO + 2) (2 772 + 77 + 1) (72 + TO + 1). 

In the 72 and to large limit this gives the rational parametrization of the singular 
variety (|E.7|) . namely (x, y) = (X(n, to), Y(n, to)), with: 

4 

72 

^(n, to) = -r— — -, Y(n, to) = X(m, n). 

(2 72 + to) 2 (2 to + 77) (77 + to) 

Of course, all these calculations can be performed with series of any finite number 
of complex variables. These (simple) calculations are only valid for series of several 
complex variables, such that the ratio of the various consecutive coefficients (see (14"4"1) ) 
are rational expressions (typically hypergeometric series). 
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6. Towards singular manifolds of Ising model D-finite system of PDEs 

One thus sees, from the previous calculations, that one can actually define, and find 
without ambiguity, the singular manifolds of D-finite systems of PDEs. The singular 
manifolds are fixed, and can (in principle) be obtained from (possibly tedious but 
well-defined) calculations from the D-finite system of PDEs. This is quite different 
from the case of generic (non-holonomic) systems of PDEs where singularities depend 
on initial boundary conditions. With the previous calculations, one can see that the 
singular manifolds can even be obtained from very simple calculations in the (selected) 
case of hypergeometric series, the singular varieties with rational parametrization being 
underlined. 

For functions of several complex variables which are not known to be solutions 
of D-finite systems of partial linear differential operators (or even partial non-linear 
differential operators but with fixed critical points), the question of defining and 
finding the singular manifolds seems hopeless. There is, however, one category 
of functions of several complex variables that emerges quite naturally in physics, 
where some hope remains, thus partially justifying, the "guessing" approach often 
performed in lattice statistical mechanics [81[2ffl[I7J[2SG3[S81[55]. These 
functions of several complex variables are the ones which can be decomposed as 
infinite sums of D-finite functions (in a typical Feynman diagram approach). The 
best example is the full susceptibility of the anisotropic square Ising model which 
has such a decomposition 37J. Let us try to find the singularity manifolds of the 
anisotropic X s : trying in a second step, to understand the singularity manifolds of 
the anisotropic full susceptibility x- 

6.1. Landau approach for the singular manifolds of the anisotropic x^ 

Finding the Fuchsian (and in fact globally nilpotent [60]) linear ODEs for the n-fold 
integrals X^' s of the decomposition of the full magnetic susceptibility of the square 
lattice Ising model is already a "tour-de-force" in the isotropic case 35 , 38 , 6TJ [521 E3] ■ 
The anisotropic X i nas a surprisingly nice factorized form (see equation (3.22) 
in [41]). It is the product of the isotropic x^ an d of a simple square- root algebraic 
function: 

((l + fcr).(fc + r)) 1/2 
X (2) (fc,r) = ^ — >—. x W(k,i), (50) 

where k — s\ S2 is the modulus of elliptic functions in the parametrization of the 
model, where the ratio r = s\/s2 is the anisotropy variable, with si = sinh2i^i, 
s 2 = sinh2if 2 , (with notations K x = E v /k B T and K 2 = E h /k B T, see (3.22) 
of [H]), and where 1) is the isotropic x^ : 

(2) J_ (l + fc 2 ).£(fc 2 )-(l-fc 2 )-if(fc 2 ) 

X { > ' 3tt' (l-fc)fl-fc 2 ) 

- i(iw' ,Fl ( [ ^ M3|i (nV)' (51) 

Beyond this surprisingly simple x^ case, obtaining a D-finite (Picard-Fuchs) 
system for x^K f° r the anisotropic square Ising model, would require too massive and 
extreme computer calculations. Furthermore, the simple "Horn calculations" detailed 
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in section ([S]) require some closed asymptotic formula (or some asymptotic formula of 
exact linear recursions) for the coefficients of the double series of the anisotropic x^ n \ 
and would require some assumption that the X s are hypergeometric series, or at 
least, that their singular part is dominated by hypergeometric series. 

However, if one is only interested in the singularities of such D-finite n-fold 
integrals, the Landau singularity approach, we have already used in the isotropic 
case, to find [3j [40] these singularities, can again, be worked out. We are not going 
to recall the details of this approach, which correspond in the anisotropic case, to 
sometimes quite tedious (algebraic) calculations. The idea, which is specific of re- 
fold integrals of some algebraic integrands, amounts to saying that the singularities 
should, in principle, be deduced only from the algebraic integrands of these integrals 
from elementary algebraic calculations [U [21 OH HOI ESS ■ 

We will display, in a following subsection (16.31) the results for the first X s after 
recalling in the next subsection a first set of fundamental singularities. 

6.2. Nickellian singular manifolds for the anisotropic x^ ' s an d zeroes of the 
partition function 

In contrast to the form factors [HI El] (M, N) , whose only singular points are 
k = 0, k = 1 and k = oo, the x W s have many further singularities. The first 
set of these singularities was found, by Nickel [121 23], to be, for the isotropic case 
(Ki = K 2 = K), located at 

cosh 2 2# - $mh2K ■ (cos(2nj/n) +cos(2Trl/n)) = 0, (52) 

with ([x] being the integer part of x): < j, I < [n/2], j = I = excluded (for n 
even, j + I = n/2 is also excluded). Equivalently ([52]) reads: 

sinh2Kj y i = Sjj, = 1/2- (cos(27rj'/n) + cos(27r//n)) 

± i/2 • [(4 - (cos(2vrj/n) + cos(27rZ/n)) 2 ] 1/2 - (53) 

These Nickel's singularities are clearly on the unit circle |s| = 1, or fc| = 1. Do 
note that this is no longer the case for the anisotropic model where Nickel's singularities 
for the anisotropic x^'s become: 

cosl^ifi • cosh 2K 2 (54) 
- (sinh2ifi- cos(2nj/n) +sinh2if 2 - cos(27r?/n)) = 0, 

with j, I = 1, 2, • ■ • , n. These (complex) algebraic curves (1541) . in the two complex 
variables s% = sinh2i^i, s 2 — sinh 2K 2 , have to be singular loci (as will be suggested 
in the following section) for the D-finite system of PDEs satisfied by the anisotropic 
(holonomic) x'"' ) ' s - 

One can rewrite these algebraic curves in k = s± ■ s 2 and r = s\/ s 2 as 

(r + fc) • (fer+1) -k-(rU±V) 2 = 0, (55) 

where U — cos(27rj'/n) and V — cos(2irl / n) . Do remark that these algebraic curves 
depend on the anisotropy variable r — s\/s 2 . We will underline this important fact 
in subsection (|6.4I) . Remarkably these curves are generically of genus-on^, not only 

f For U = V (as well as U = —V, U = ±1, V = ±1) the curves are genus-zero. For instance, 
for U = V, they read (r ± l) 2 k ■ U 2 - (r + k) ■ (k r + 1) = 0. 
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when U — cos(27rj/n) and V = cos(27rZ/n), but for any fixed value of U and V. 
Their j-invariant [101 [HS] readijj 

256 (^ + ^-^-^-^ + 1)3 

(I/ 2 - l) 2 (U 2 - l) 2 (CI 2 - V 2 ) 2 ' 1 ^ 

We thus see that we do have a two-parameters family of elliptic curves. 

These elliptic (or rational) curves (|54|) accumulate with increasing values of n, in 
the same way Nickel's singularities (f52|) accumulate on the unit circle \s\ = 1, in a 
certain (real) submanifold S of the two complex variables s±, S2 (four real variables). 
However, this "singularity manifold" S is not a codimcnsion-one (real) submanifold 
(like the unit circle |s| = 1 in the s-complex plane), but actually a codimension zero 
submanifold, as can also be seen on various analysis of complex temperature zeroes 
(see for instance [69J [70l [TTJ [72]) ■ Note that this "singularity manifold" becomes very 
"slim" near the (critical) algebraic curve k = si S2 = 1 (see for instance the region 
near the real axis of figures 1, 2 and 3 in [69] ). 

In the isotropic case, we actually obtained [35l HH [39l H2 [63l [73] the linear ODEs 
satisfied by the first x^' s j for n = 3, 4, 5, 6 and, thus, of course, the corresponding 
ODE singularities. Furthermore, we also performed a Landau singularity approach 
that enabled us to obtain, and describe, the singularities for all [40] the X s - 
These exact results show, very clearly, that there are (non-Nickellian) singularities 
inside the unit circle and outside the unit circle (see Figure 1,2,3 and 4 in [40)). On the 
figures of [40] it is easy to get convinced that the accumulation of these non-Nickellian 
singularities will probably be a dense set of points inside the unit circle and (by 
Kramers- Wannier duality) outside the unit circle. These non-Nickellian singularities 
are given in terms of Chebyshev polynomials of the first and second kind (see equations 
(28) and (29) in [40). Upgrading these slightly involved exact (Chebyshev) non- 
Nickellian results [40] for the isotropic model to the anisotropic model is, at the present 
moment, probably too ambitious. 

Let us simply try, using the previous Landau singularity approach, to provide, may 
be not an exhaustive description of all the singularities for the anisotropic case, but at 
least, the exact expression of all the singular manifolds (Nickellian or non-Nickellian) 
for the first anisotropic x s - 

6. 3. Singular manifolds for the first anisotropic x^ 

The Landau singularity approach detailed in [3l [40] for the isotropic x^' s 01 the 
square Ising model, can easily be generalized to the anisotropic x^' s - We are not 
going to explain here the details of these (slightly tedious) calculations which arc 
basically the same as in [21 HO] mutatis mutandis. The calculations being slightly 
involved we just give the results for the first X^' s - 

The singularities of x^ an d X^ read respectively in k and r: 
Sing(x {3) ) = (k 2 - 1) • (3 kr + r + 4 k 2 ) ■ (k 2 r + 3 kr + 4) • (k 2 r + r + k) 

x (3r 2 k-r-k-k 2 r) ■ (4 + 3 kr + 4 k + 4 k 2 ) (r + k) {kr + 1), (57) 
Sing{x {4) ) = (k 2 - 1) ■ (kr + 1 + k 2 ) ■ (3 r 2 k - r - k - k 2 r). (58) 

\ This rational expression H56H of U and V is nothing but relation (36) in I66| with J x /Jz = U, 
Jy/Jz = V. This rational expression remarkably factorizes for many Heegner numbers 67. 68 
(complex multiplication cases): j = 12 3 , 20 3 , (— 15) 3 , 2 • 30 3 , 66 3 and selected quadratic values 
of j-invariant, like j 2 + 191025 j - 495 3 = or j 2 - 1264000 j - 880 3 = 0. This (partially) 
explains the occurrence in (1">41i of several complex multiplication cases (for instance U = cos(27r2/8), 
V = cos(2vr/8) which give j = 1728). 
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In order to compare these results with our previous exact results for the isotropic 
model, which were given [3SJ |3H [02] in the (quite natural for such n-fold integrals) 
variable [38, 62 w, let us rewrite these results in r and w = s/(l + s 2 )/2, where, 
now, s — (si S2) 1 / 2 : 

Sing(x^) = (w 2 - 1) • w 2 ■ (r 2 - 4r + 4 + 3w 2 r 2 - Aw 2 r + 16 w 4 r) 2 
x (1 + 4 w 2 r - 2 r) 2 (3 r 2 - 1 - 4 w 2 r + 2 r) 2 

x (3r-4+16w 2 ) 2 -(l + 4w 2 r-2r + r 2 ) 2 , (59) 
Sing{x {4) ) = w 2 ■ (w 2 - 1) • (4w 2 - 2 + r) 2 • (3r 2 - 1 - Aw 2 r + 2r) 2 . (60) 

Note that the complex multiplication points of the isotropic case [30], namely the 
roots ofl + 3fc + 4fc 2 = and k 2 + 3 k + 4 = 0, come from the Sing(x^ 3 ') factor 

r 2 -4r+4+3 w 2 r 2 - 4w 2 r + 16 w 4 r, (61) 

in (|59p. or equivalently with (fc, r), the two factors in ([57]) : 

(3fcr +r +4fc 2 ) • (fc 2 r + 3 fcr +4), (62) 

The vanishing condition of (I61[) corresponds to a genus-zero curve, its rational 
parametrization being: 

" 2 + l -4 

w = ~~ « 1 r = ^5 — 7 — ^ TT. (63) 

2u u 2 -(it 2 + 3) y J 

Note that Sing(x^') and Sing(x^) have a non-trivial gcd (respectively in fc, 
then w): 

gcd(Sing( X {3} ), Sing{ X {i) )) = (k 2 - 1) ■ (3 r 2 k - r - k - fc 2 r), 
gcd(Sing( X {3) ), Sing( X (4) )) = w 2 ■ (1 - w) (1 + w) ■ (3 r 2 — 1 — 4 iu 2 r + 2 r) 2 , 

the last algebraic curve 3r 2 k — r — k — k 2 r = 0, is a genus-one curve. A way to 
understand, in the anisotropic case, the emergence of singular algebraic curves shared 
by several X ^ 's ( n even and n odd) amounts to noticing that these curves actually 
reduce, in the isotropic limit, to k — 1, the singular variety of the partition function 
of the anisotropic model. 

The fact that the singular curve 3r 2 k — r — k — k 2 r = 0, together with the 
Nickellian algebraic curves (|54j) . (|55j) . are not genus-zero (as all the genus-zero curves 
of section (g]), like ([2"Qjl. ([3~4"]). as well as the ones displayed in |Appendix E[ see (|E.3p . 
l|E.7p ), show that the series for the anisotropic X s cannot be hypergeometric series 
in the variables k and r (see section ([SJ). 

It would be interesting, before trying to generalize the Chebyshev polynomial 
formula [10J for the non-Nickellian singularities of the isotropic model, to the 
anisotropic one, to accumulate, with this Landau singularity approach, more non- 
Nickellian algebraic curves in the anisotropic case. Recalling the systematic emergence 
of elliptic curves (see ([55]) ) for the Nickellian algebraic curves, it would be interesting 
to systematically look at the genus of these singular curves, to see if higher genus 
curves are also discarded for the non-Nickellian algebraic curves. 

It would be also interesting to confirm these Landau singularity calculations, with 
differential algebra calculations. Even with the last progress performed by Koutschan 
on the creative telescopic method [74] [75], getting the (Picard-Fuchs) system of PDEs 
satisfied by the several complex variables series of the anisotropic x^'s corresponds, 
at the present moment, to too large calculations (even for the anisotropic x^)- 



Singularities of lattice models 



20 



However, if one considers particular anisotropic subcases ( s 2 = 3si, s 2 — 5s^, ...), 
obtaining the corresponding ODEs for the anisotropic x*- 3 -*, in the unique complex 
variable, could be imagined using the creative telescopic method [7H [75], or even, 
from series expansion as we did in the isotropic case |35j . 

6.4- Singular manifolds and the anisotropy variable 

For experts of Yang-Baxter integrability, the fact that the singularities varieties, 
namely the Nickellian elliptic curves (f55|) . or the non-Nickellian rational curves l[62jl . 
do depend on the anisotropy of the model may come as a surprise. Indeed, within the 
Yang-Baxter integrable framework, and as a consequence of the existence of families 
of commuting transfer matrices (row-to- row, diagonal or corner transfer matrices), 
one used to have many quantities like the order parameter, the eigenvectors of row- 
to-row or corner transfer matrices, which are independent of the so-called "spectral 
parameter" (the parameter that enables to move along each elliptic curve). The 
selected quantities depend only on the modulus k of the elliptic functions. Along 
this line, one certainly expects the singular manifolds, which are highly symmetric, 
"invariant" and "universal" manifolds 1151 116) . to be also independent of the 
spectral variables. With the previous variables k and r, the singular manifolds should 
just depend on the modulus fc, and not on the anisotropy variable r (related to the 
spectral parameter). The surprise is that the singular manifolds do depend also on 
the anisotropy variable r, and thus on the spectral variable. 

The X s are known |64) to be an infinite sum of form factors C^ n \N, M): 

x (n) = £^C<(n) (iV)M)) (64) 
M N 

this relation being inherited from the fact that the full susceptibility is the sum of all 
the two-point correlation functions |64j . 

Recalling the simplest (nearest neighbour) correlation function C(0, 1), it 
reads [76] in the anisotropic cas 

C(0,1) = — • (^r) • ((1 + *r)-n(-Ar, A) - #(*)), 

where, again, k = s\ s 2 is the modulus of the elliptic functions parametrizing the 
model, and r is the ratio r = s\/s 2 and where Tl(x, y) is the complete elliptic 
integral of the third kind. 

The singular manifolds correspond to the singular points of the complete elliptic 
integrals of the first and third kind, namely k = 0, k = 1 and k — oo. Therefore 
they depend only on the modulus k in the elliptic parametrization of the model. 

The form factors have been seen to be solutions of linear differential equations 
associated with elliptic functions |641 165) . Consequently, their singular points 
correspond to the singular points of the complete elliptic integrals of the first or 
second kind E or K, namely k = 0, A; = 1 and k = oo. The generalization to the 
anisotropic case has been sketched in [76]. One expects the results to be polynomial 
expressions of the complete elliptic integrals of the first (or second) and third kind, 
yielding again, singular manifolds which depend only on the modulus k, and are 
actually k = 0, k = 1, or k = oo. 

Finite sums of correlation functions or form factors, certainly have k = 0, k = 1 
or k — oo as singularities, even for the anisotropic model. However, the anisotropic 

% We use the maple notations for IT and K. 
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^,(n) > g are sums f arL infinite number of form factors. One cannot try to deduce 
the singular points of these infinite sums x'™^ s from the singular points of the form 
factors. The X s are ' m f ac ^: quite involved "composite" quantities with no simple 
combinatorics interpretation (like being the sum over graphs of a certain type). It 
is worth noting that exploring all the algebraic singular curves for all the X s : 
condition k — 1 always occurs for all the X s - 

The previous results provide a quite interesting insight on the "true mathematical 
and physical" nature of the x^' s: they are quite involved "composite" quantities, 
their singularities being drastically different from the ones of the C^- n \N, M) form 
factors [64 1 165 ] , 

In the isotropic case, strong evidence has been given [39] [40] [41] [42] [43] that 
the full susceptibility x nas a natural boundary corresponding to the accumulation of 
singular points on the \k\ — 1 unit circle, thus discarding a common wisdom that "of 
course" the singularities of the partition function are the same as the singularities of 
the full susceptibility. 

By analogy with the situation encountered in the isotropic case, we are going to 
have an accumulation of singular curves densifying the whole parameter space (two 
complex variables s\ and S2, i.e. four real variables). The equivalent of the unit circle 
is now, a codimension-zero manifold in the four real variables parameter space, which 
disentangles two codimension-zero domains in the parameter space. Is it the singular 
locus for the full anisotropic susceptibility x ? Do we have here a generalization of 
the concept of natural boundary for several complex variables ? If the answer to 
the question of the location of the singularities of non-holonomic functions seems to 
be dependent of the decomposition of the non-holonomic function in infinite sums of 
holonomic functions, is it simply well-defined ? 

All we can reasonably say is that, probably, and in the same way as in the isotropic 
case, the double series for the x^' s are n °t singular in one domain (the equivalent 
of the inside of the unit circle), and one probably has the same result for the full 
anisotropic susceptibility x- 

6.5. Anisotropic models: n-fold integrals of several complex variables 

In the anisotropic case, the X^' s are "-fold integrals of several complex variables. 
After Kashiwara and Kawai [28] , we do know that these "functions" of several complex 
variables are holonomic. Let us restrict to the case, we often encounter in physics, 
where the integrand is an algebraic function of these several complex variables (and 
of the integration variables). In contrast with the one complex variable case, the 
holonomic character, here, corresponds to an extremely rich structure: the solutions 
of the over- determined system of linear PDEs correspond to a finite set of solutions (for 
one complex variable this is obvious), and the singularities, which are no longer points 
but manifolds, are fixed algebraic varieties (for one complex variable this is obvious). 
Furthermore these operators are globally nilpotent (the holonomic functions can, in 
this "Derived from Geometry" framework [73 023 , be interpreted as "Periods" of an 
algebraic variety closely related to the integrand). We have many other remarkable 
properties. For instance, the operators are often (always ?) homomorphic to their 
formal adjoint (this is related to the occurrence of selected differential Galois groups). 
All these remarkable properties correspond to a differential algebra description of 
these structures. Finally, we have also other properties of more arithmetic and 
algebraic geometry nature. The series expansions of these holonomic functions are 
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often globally bounded [50 , which means that they can be recast (after rescaling) into 
series expansions of several variables with integer coefficients. This raises the question 
of the '' 'modularity" in these problems [77] [75]. Along this "modularity" line, beyond 
the occurrence of many modular forms [791 180] , we also see the emergence of Calabi- 
Yau ODEs. From a differential algebra perspective, the emergence of Calabi-Yau 
structures [81] is not clear. In some integrability framework, the argument that Calabi- 
Yau manifolds are, after K3 surfaces, the "next" generalization of elliptic curves, 
remains an insufficient and much too general argument. 

Let us inject, beyond the differential algebra description of these structures, some 
birational algebraic geometry ideas. In lattice statistical mechanics, the models defined 
by local Boltzmann weights depending on several complex variables, are known to 
have, generically, an infinite set of birational symmetries generated by the combination 
of the so-called inversion relations 82, S3J. 

It has been shown that n-fold integrals like the X^' s of the Ising model present 
some nice inversion relation functional equations in the anisotropic case [84j (several 
complex variables): 



inherited from the same inversion relation functional equation on the full anisotropic 
susceptibility. 

Since the previous ideas underline the crucial role of the integrand of the n- 
fold integrals as the algebraic variety from which "everything" , in principle, can be 
deduced HH1 [75] Hi], it is interesting to see if this integrand, itself, is not going 
to be invariant (resp. covariant) by these birational involutions (and, thus, by the 
composition of these birational involutions) when we keep the integration variables 
fixed. One can verify that this is actually the case for the integrand of the anisotropic 
)(W's of the Ising model. 

Unfortunately, the group of birational transformations of the Ising model is a 
finite set of transformations. However, for generic models, one can easily imagine 
to be in a situation where the integrand of the n-fold integrals of several complex 
variables emerging in these models, will be invariant (resp. covariant) by an infinite 
set of birational transformations [S] . 

We will thus have a natural emergence (in lattice statistical mechanics) of 
algebraic varieties with an infinite set of birational symmetries [5]. These algebraic 
varieties have zero canonical class, Kodaira dimension zero. We, now, understand the 
emergence of Calabi-Yau manifolds in these problems: Abelian varieties and Calabi- 
Yau manifolds (in dimension one, elliptic curves; in dimension two, complex tori and 
K3 surfaces) have Kodaira dimension zer<U 

One can expect that the singular varieties (like ([9]) or (|12p ) will have to be 
invariant by the (generically infinite) set of birational transformations generated by the 
inversion relations. When the singular manifolds are algebraic curves, the existence of 
a (generically infinite) set of birational automorphisms for the algebraic curves implies 
that the curves are, necessarily, genus zero or one [5] . This enables to understand!^ the 
emergence of remarkable structures like the two-parameters family of elliptic curves 

f Zero canonical class, corresponding to admitting flat metrics and Ricci flat metrics, respectively. 
§ Cum grano salis: in the (free-fermion) Ising case the birational transformations generated by the 
two inversion relations form a finite set 1901 191| . which allows, in principle higher genus curves. One 
must imagine the Ising model as a subcase of a larger model with n-fold integrals, where one would 
recover a (generic) infinite set of birational transformations. 



X 




(65) 
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(|55p. Actually this is the way many singular varieties have been discovered on many 
lattice statistical mechanics models (see [TH [THl HH H] ) • This birational invariance 
fits quite well with the interpretation of the singular variety (|12l) , as the discriminant 
of a two-parameters family of Calabi-Yau 3-folds. 

7. Conclusion 

In the theory of critical phenomena (renormalization group, etc), singularities are often 
seen as fixed points of a "dynamical system" called renormalization [S6J , and one takes 
for granted, with a (lex parsimoniae) simplicity prejudice, that these singularities are 
isolated points, or smooth manifolds (hopefully algebraic varieties [H|9l[Tp ( 87, 88,(89] if 
one has an integrability prejudice as well). In the theory of discrete dynamical systems, 
a totally opposite prejudice exists like the belief in a frequent occurrence of strange 
attractors for the set of fixed points of many "dynamical systems" . Singularity theory 
in mathematics, and in particular Arnolds's viewpoint |27) . are a perfect illustration 
that the set of singular points should actually correspond to much more involved 
manifolds than what is expected in the mainstream doxa of critical phenomena. 

We have performed some kind of "deconstruction"0 of the concept of singularities 
in lattice statistical mechanics. The sets of singularities are much more complex sets 
of points than what physicists tend to believe (see Figures 1, 2, 3, 4 of [4"0]h 

The mathematician's viewpoint that singularities are much more complex than 
what physicists believe with their (lex parsimoniae) simplicity optimism, is the correct 
viewpoint. On the other side, the mathematician's viewpoint that nothing serious 
and/or rigorous can be done with several complex variables is too pessimistic: within 
that viewpoint, singularities are seen as too involved to analyze, impossible to localize 
(of course outside the hypergeometric series framework), or simply, a not well-defined 
concept. Even in the case of several complex variables, many singular manifolds 
conjectured by physicists, in particular F.Y. Wu E] [10], turned out to be true 
singular varieties of lattice models, because physicists are (sometimes without being 
fully conscious) often working with holonomic (D-finite) functions of several complex 
variables. 

Focusing on the full susceptibility x °f the (anisotropic) Ising model and on the 
holonomic X s j we have obtained singular manifolds of the linear partial differential 
systems of the X s - The fact that these singular manifolds do depend on the spectral 
parameter of this Yang-Baxter integrable model is a strong indication that these x'™" 1 ' s 
are highly composite objects (even if the exact expression of these singular varieties 
remains simple enough for the first X ' s )' Furthermore, the fact that most of these 
singular manifolds are not genus-zero curves show that the series of the anisotropic 
t^W's, despite all their remarkable properties, cannot be reduced to hypergeometric 
series. 

In the case of the full susceptibility x °f the (anisotropic) Ising model, we seem to 
have the following situation: among the quite large, and rich, set of singular varieties 
of the linear ODEs of the x^' s i there is a restricted set (see (|54"]) . ((55]) ) of singular 
varieties which actually corresponds to zeroes of the (anisotropic) partition function, 
and, in the same time, corresponds to singularities of the linear PDEs of the X s - 
This set could correspond (by analogy with the isotropic case) to singularities of the 
series expansions of the x^' s - A fundamental idea to keep in mind is that it is 

X Using Derrida's wording. 
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crucial to make a difference between the singularities of the (series expansions of the) 
D-finite functions, and the singularities^] of the linear partial differential systems for 
these functions. 

It would be interesting to see if, inside some reasonable theoretical physics 
framework, similar result^ can also be obtained for other non-holonomic functions of 
several complex variables that decompose into an infinite set of holonomic (D-finite) 
functions. 
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Appendix A. The nine formal solutions of the Picard-Fuchs system 
"above" the Calabi-Yau ODE © 

Let us find the "formal solutions" around (x, y) — (0, 0), of the PDE system ([5]) 
"above" the Calabi-Yau ODE Q. One plugs, in ©, the series 

U 3,k{x,y) ■ Hxf \n( y y- k , (A.l) 

j=0 k=0 

where T-Lj t k(x, y) are series in x and y and solves the system term by term. Collecting 
on the non fixed coefficients, one finds So = Ho(x,y) and 

Si = H {x, y) ■ ln(x) + Hi{x, y), S 2 = H (x, y) ■ Info) + H 1 (y, x), 

S 3 = H (x, y) ■ ln(x) 2 +2H 1 (x, y) ■ ln(a:) + H 2 (x, y), 
Si = H (x, y) ■ Info) 2 +2H 1 (y, x) • Info) + H 2 (y, x), 

5 5 = H Q (x, y) ■ \n(x) ■ Info) + Hi(y, x) ■ \n(x) + Hi(x, y) ■ Info) + H 3 (x, y), 

5 6 = H (x, y) ■ ln(x) 2 • Info) + 2H 1 (x, y) ■ \n{x) ■ Info) + H x {y, x) ■ \n{x) 2 

+ 2H 3 {x, y) ■ hx{x) + H 2 (x, y) ■ Info) + H 4 {x, y), 

5 7 = H (x, y) ■ ln(ar) • ln(y) 2 + 2H x {y, x) ■ ln(x) • Info) + H x {x, y) ■ ln(y) 2 
+ 2H 3 (x, y) ■ Info) + H 2 (y, x) ■ ln(x) + H 4 (y, x), 

5 8 = H Q {x, y) ■ ln(x) 2 • ln(y) 2 + 2H 1 (y, x) ■ ln(x) 2 • Info) + 2H 1 {x, y) ■ hi{x) ■ ln(y) 2 

+ 4i7 3 (x, y) ■ \n{x) ■ ln(y) + H 2 (y, x) ■ ln(x) 2 + H 2 (x, y) ■ ln(y) 2 

+ 2H 4 (y, x) ■ \n{x) + 2H 4 (x, y) ■ Info) + H 5 (x, y), (A.2) 

^ The singular manifolds seem to have, in the case of n-fold integrals of algebraic integrand, a 
projective invariant interpretation as discriminant of the algebraic varieties associated with the 
integrand. 

f With the problem that the results seem, at first sight, to depend on the decomposition in an infinite 
sum of holonomic functions. 
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where (only the first terms of the series are given) 

H (x,y) = 1 + 6(x + y) + (90 (x 2 + y 2 ) + 720 xy) + 

/ 5 13 1323 
H x {x,y) = (15z+33y) + (— x 2 + 3132 xy + — y 1 

/279 

H 2 (x,y) = (108 y -18a;) - ( — x 2 -%120xy - 3654y 2 

/2709 2709 
H 3 {x,y) = 9- (x + y) + ( —x 2 + 3960 xy + — y' 

^x 2 +11178xy +—y 2 
( 14931 

H 5 (x,y) = 324- (x +y) - ( ~^{x 2 + y 2 ) - 6912xy 

There are nine solutions for the system ([5]). One notes that Ho, H 3 and H§ 
are symmetric in x, y, while Hi, i?2 and Hi are not symmetric in x, y. For the 
formal solutions, So, S5 and Ss are symmetric in x, y, and the six others are pairwisc 
symmetric. These nine independent formal solutions are solutions of the PDE system 
(|5"T). and thus of the order-nine differential operator VL X and its (x, y)-symmetric £l y . 

Note however, that the linear differential operator VL X has been constructed from 
the PDE system ([5]) and factorizes as written in (fl7|) . it, then, might be that Ho(x,y) 
is a solution of only the right factor operator Lq(x, y). Indeed, plugging a series 

into -Lg(a;, y) and solving term by term, one obtains (up to the overall co,o), the double 
hypergeometric series Hq(x, y). The solutions of L§(x,y) can be expressed in terms 
of the previous formal solutions (| A.2|) : 

S 4 

So, Si, S2, S3 — S4, 55 H — — , ^6 + Sj. (A. 4) 

Appendix B. Factorization (|1T|> of the order-nine operator Cl x 

The order- nine operator £l x of subsection fl3.2[) factorizes (see (TTTl) ) into three order- 
one operators and the order-six operator L 6 (x, y): 

1 



L 6 (x, y) = — • ^ p n (x, y) ■ D n x , (B.l) 

P6[X > y > n=0 

The three order-one operators are encoded by three rational functions of x and 
y, namely ri(x, y), r2[x, y) and f 3 (x, y). These polynomials factorize (see (|T9l) ) 
and thus the fi(x, y)'s reduce to the expressions of four polynomials with integer 
coefficients Vg(x, y), V§(x, y), q\ and q 2 , where Vg(x, y) is the polynomial of the 
apparent singularities of the order-nine operator il x , and where V&(x, y) is the 
polynomial of the apparent singularities of the order-six operator Lq(x, y). 
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These polynomials read: 
T 9 (x,y) = 2 4 -3 18 -x 6 -2-3 16 • (31951 + 1602072 y) ■ x 5 
+ 3 13 • (14397329 + 913784868 y+ 17712588816 y 2 ) ■ x 4 
+ 3 9 • (2986814425 + 60616383939 y- 1350750590172 y 2 

- 24695209500192 y 3 ) ■ x 3 

+ 3 7 • (5310925151 - 333452529387 y- 14254789072275 y 2 

+ 241096254564492 y 3 + 7702353325801296 y 4 ) ■ x 2 

- 81 • (27 y- 1) ■ (39319888296092688 y 4 + 122020942792986 y 3 

- 111685613173821 y 2 + 22118310900 y + 86524357339) • x 

+ 2 4 -5 3 • (10827 y + 364) 3 - (27y-l) 3 , (B.2) 

V 6 {x, y) = 387420489 • (x 2 - 142 xy + 343 y 2 ) ■ (x + y) 4 

- 43046721 • (x + y) • (89 x 4 - 196 y 4 - 823 xy 3 + 13287 x 2 y 2 - 3493 x 3 y) 
+ 1594323 • (3482 x 4 + 662 xy 3 + 2972 x 3 y - 427 y 4 + 25365 x 2 y 2 ) 

+ 19683 • (33307 a; 3 - 1784 y 3 - 14487 xy 2 + 44904 x 2 y) (B.3) 

- 2187- (27394 x 2 ~ 88 xy- 671 y 2 ) + 162 • (1325 x + 242 y) -1331, 

<7i = 4 • 3 18 • (x 6 + 113061462 xy 5 + 4560 x 5 y - 8876482 x 3 y 3 + 284847 x 4 y 2 
- 52726107 xV + 28140175 y 6 ) 
+ 3 16 • (4108 x 5 - 11112875 x 3 y 2 + 587276 x 4 y - 105291883 xy 4 

+ 14516200 y 5 - 4491914 x 2 y 3 ) 
+ 3 13 - (198311a; 4 - 370624786 xy 3 + 6765614a; 3 y 

- 130714000y 4 + 116112144xV) 

+ 3 9 • (18879841 x 3 - 64727000 y 3 + 773936148 xy 2 + 17519934 x 2 y) 

- 3 7 • (45403057 x 2 - 221205178 xy - 141045500 y 2 ) 

- 567- (22002263 x- 1112800 7/) - 145745600, 

q 2 = 774840978- (x 6 - 841926 xy 5 - 462x 5 y - 341728 x 3 y 3 + 32721 x 4 y 2 
+ 810681 x 2 y 4 + 98245 y 6 ) 

- 43046721 ■ (223 x 5 + 54121 x 3 y 2 - 47245 x 4 y - 613336 xy 4 

- 68810 y 5 + 20707 x 2 y 3 ) 

+ 1594323- (22489 x 4 + 1358236 xy 3 + 304861 x 3 y 

- 250820 y 4 - 1645923 x 2 y 2 ) 

+ 19683 (415049 x 3 - 505660 y 3 - 4725138 xy 2 + 65103 x 2 y) 
+ 10935- (229157 x 2 - 163880 xy + 68006 y 2 ) 
+ 162 • (492079 x + 45925 y) -440440. 
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Appendix C. Alternative linear differential operator for the double 
hypergeometric series 

Recalling the double hypergeometric series ([lj, Hq(x, cx) is solution of an order-six 
c-dependent linear differential operator 

W 6 = (1 +162-(c+l)-x) x 

(f - 81 • (c + f ) • x + 2187 • (c 2 - 7c + 1) • x 2 - 19683 ■ (c + I) 3 ■ x 3 ) ■ x 4 ■ D% 

+ •■• (C.l) 

In the c = 1 limit, this order-six operator becomes the direct sum of the order- 
two linear differential operator 

2 - 3a;- (30 + 1) • (30 + 2), 

with the hypergeometric function solution 

2 *i([i 2 -]M-27x), (C.2) 

and of the order- four Calabi-Yau ODE ([JJ , with the analytic solution ([3]) , which can 
be written as the Hadamard product [9"2"] : 

•*([!, -27*) * • §], [l]; -j^)). 

In the (less natural) c = limit, this order-six linear differential operator is the 
product of homomorphic operators: 

W 6 {c = 0) = N 2 -M 2 -L 2 , (C.3) 

where L 2 has the hypergeometric function solution 

2 F x f[i §],[!]; 27 A (C.4) 



L 3' 3 

In the c — y co limit, this order-six operator degenerates into the direct sum: 
(3 • + 1) © (3 • 9 + 2) © (3 ■ 9 + 4) © (3 • 6 + 5) © (3 • + 7) ® (3 • + 8). 

Appendix D. Another series of two complex variables 
Appendix D.l. Double hypergeometric series 

Without the factor 64, the results for (f3"U| in subsection (|4.2[) correspond to the double 
hypergeometric series 

00 00 < \3 (a\3 



{a) 3 n ■ {fif m ■ 09') 



n „,m 



where (a)„ is the usual Pochhammer symbol. The double hypergeometric series 
K(x, y) is a Kampe-de-Feriet function [4^14911461147] 

^o.o ([£'], [«, «, a], [/3, /?, /?]; [7, 7, 7], - -; x, y). (D.l) 
The singularity varieties of (|D.1I) arc independent of the parameters a, /3, /?', 7, 
and are x ■ (1 — x) ■ (1 — y) ■ (y — x) = 0, together with 

y 2 x 2 -2xy{y + x) + (x - y) 2 = 0, (D.2) 

in agreement, in the a = j3 = f3' = 1/2, 7 = 1 limit, with (|34|) . taking into account 
the rescaling (x, y) — > (64 x, 64 y). 
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Appendix D.2. Other double hypergeometric series 
Introducing the other double hypergeometric series 

K2(x,y) = EE {a) \'^'fm^ ' xn ' ym - (D ' 3) 

It is also a Kampe-de-Feriet function [46j |47l HH1 149] 

C,o M (^']. [a, • • • , a], • • • , ^]; [7, • • • , 7], - "5 2/)- (D.4) 
Let us restrict, in the following, to a = /3 = ft = 1/2 and 7 = 1. 

The singularity varieties of the PDE system are actually different from (|D.2|) and 
depend on M. For M = 2 and M = 4, they read respectively: 

(x + y) 2 - x 2 y 2 = 0, {x + y- xyf + 27 x 2 y 2 = 0. (D.5) 

More generally, for M an even integer, besides the conditions x- (1 — x) ■ (1 — y) = 0, 
the singular manifold reads an algebraic curve of parametrization 

/ —t \ M-l 

x = t*-\ y = (y^) , (D.6) 
or equivalently 

" (a+f) , » = (a-f) , (D.7) 

that can be thought as a "Fermat-like" curve: 

a;T^w + yriu = 1, (D-8) 

For M = 3, we have (|D.2|) and for M = 5, we have (besides the conditions 
x ■ (1 — x) • (1 — y) ■ (y — x) = 0) the singular variety 

(x + y + x j/) 4 — 136 x 2 y 2 ■ (x + y + xy) — 8 x y ■ (x + 1 + y) (x 2 + y 2 ) 

-8x 2 y 2 ■ (x + y)(xy- 1) = 0. (D.9) 

More generally, for M an odd integer, besides the conditions x- (1— x)- (1— y)- (y—x) — 
0, the singular manifold reads an algebraic curve of parametrization 

x = t M -\ y = (^-L) , (D.10) 

or equivalently 

x = -o + v ) . 2/ = -o -«) . (D-ll) 



2 / V 2 

that can be thought as a "Fermat-like" curve: 

X T=M + yl^M +1 = 0. (D.12) 

Appendix D.3. Differential operators restricted to singular varieties 

Let us restrict to the singular variety (|D.2|) for M = 3, using the rational 
parametrization (|D.10[) . that is (x, y) — (t 2 , (t/(l — t)) 2 )- The double series 
expansion (|D.3[) becomes a series expansion in the t variable which is solution of 
the order-four linear differential operator (D t — d/dt): 

Ci = t 3 - (t-l)(2t+l)(t + 2)(t 2 + t + l) 2 (t + lf- 

+ 2t 2 ■ (t 2 + t + l) ■ (t+ l) 3 ■ cs(t) ■ D 3 t +t ■ (t+ l) 2 • c 2 (t) ■ D 2 t 

+ 2 (* + !)• ci(i)- A +2t • (t + 2)(i 2 +t+l) 4 , (D.13) 
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where 

c 3 (t) = 10i 6 + 32i 5 + 39< 4 + 20i 3 -m 2 -24i-6, 

c 2 (t) = 50i 9 + 243i 8 + 588i 7 + 903i 6 + 885i 5 + 501i 4 + 33i 3 -174i 2 -99i-14, 

ci(t) = 15 i 10 + 82 i 9 + 228 + 411 1 7 + 531 i 6 + 513 1 5 + 333 

+ 99f 3 - 12 1 2 - 12 1 - 1, 

This "critical" order-four operator C4 is such that its exterior square is a linear 
differential operator of order five (and not six as it should be for a generic order-four 
operator) . This condition that the exterior square is of order five is called the "Calabi- 
Yau condition" : it is one of the conditions defining Calabi-Yau ODEs [44j EH [SU [95] . 
Related to this exterior square condition one also has the property that this order- 
four operator C4 is homomorphic to its adjoint, up to a conjugaison by the polynomial 
(x + l) 3 (x 2 +X + if. 

Note that the limit y = x, yielding to the Calabi-Yau operator (|28p (also such 
that its exterior square is a linear differential operator of order five), is actually a 
singular limit of the Picard-Fuchs system. 

Similarly, let us restrict to the singular variety (|D.5|) for M — 2, using the 
rational parametrization (ID.6|) . namely (x, y) = (t, —t/(l — t)). The double series 
expansion (|D.3[) becomes a series expansion in the t variable which is solution of the 
order-three linear differential operator (Dt — d/dt): 

C = Df , 3. (3t-22. l.m'-iet-M , 1 1-2 



* 2 t(t-l) ' 4 {t-l) 2 -t 2 * 8 t-(t-l) 3 ' 
This "critical" order-three operator C3 is such that its symmetric square is a linear 
differential operator of order five (and not six as it should be for a generic order-three 
operator) . Related to this last property one also has the property that this order-three 
operator C3 is homomorphic to its adjoint, up to a conjugaison by the rational function 
l/x 2 /(x- 1). 

This order-three operator C3 is, in fact, exactly the symmetric square of 

16 1 • (t — l) 2 • D 2 +8 • (3t- 2) ■ (t- 1) ■ D t +t, (D.14) 

which has (1 — i) 1 / 4 • Ki^t 1 ! 2 ) as a solution (K is the complete elliptic integral of the 
first kind). 

Let us now restrict to the singular variety (|D.5|) for M — 4, using the 
(alternative) rational parametrization 

x = 8t, y = -j-^gf- (D.15) 

With this parametrization the double series expansion (|D.3I) becomes a series 
expansion in the t variable with integer coefficients. It is solution of an order-eight 
operator, its symmetric square is of order 35 (and not 36 as it should be genericalljij]) . 
For M — 4 the double series can also be resummed in one variable and rewritten 

as 

E 45^,10 • ^([~, ~ 5.5.™+ \], [m + 1, m + 1, m + 1, m + l];x) • y m 

m = 

corresponding to the identity 

(2m)! 5 /(l/2) 4 .(m + l/2)„x (1/2) 4 ■ (l/2) 4 m ■ (l/2) m+n 



4 5m • to! 10 V n\-(m + l)i J n! • m! • (l) 4 l+Il 

f Its exterior square is order 28 as it should for a generic order-eight operator. 
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More generally one has the identities 

(2m)! M (l/2)£ 



(D.16) 



4Mm . ™|2M m! M ' 

and 

(2my M+1 / (l/2)y-(m + l/2) n x = (1/2)^ • (l/2)ff • (l/2) m+n 

4 (M+l)m . m !2(M+l) V n !.( m + l)M ^ n | . m ] . (l)M +n 

and the alternative writing of the the double series (ID.3|) . for a — (3 = (3' = 1/2 
and 7 = 1, as 

^ (2m)! M+1 



L^i 4(M+l)m . m |2(M+l) 
m = 



(D.17) 



i^m([^, ■•• , m + i], [m + 1, •■• , m + l];x) ■ y" 



in- 



Let us now restrict to the singular variety y — x. For M = 4 and M = 5, 
the double series expansion (|D.3[) becomes a series expansion in a; that is solution of 
an order-six linear differential operator. For M — 4, this order-six operator is such 
that its symmetric square is of order 20 (instead of the order 21 one could expect 
generically). For M = 5, this order-six operator is such that its exterior square is of 
order 14 (instead of the order 15 one could expect generically). 

Appendix E. More Picard-Fuchs systems above Calabi-Yau ODEs 

Appendix E.l. More Picard-Fuchs system with two variables 

Another example is the two-variables Picard-Fuchs system "above" the order-four 
Calabi-Yau operator (see ODE number 18 in [44] ) ) 

e 4 - 4 x • (3 e 2 + 3 e + 1) • (2 e + i) 2 

-4j: 2 • (46> + 5) • (46> + 6) • (40 + 2) • (40 + 3) (E.l) 
= (1 -64a;) ■ (1 +16 a;) ■ x 4 ■ D* + ■■■ 
The Picard-Fuchs system corresponds to the double series [36] 



V V 



(n + to)H (2m + 2n)\ 



^ ^ n! 4 ml 4 

n— m—0 



x y 



(2m)! /. 1. . . 

= 2^ ^-ff ■ *Fs([m +l,m + l,m +1, m +-], [1, 1, 1]; Ay 

m=0 

= 1 + 2 • (x + y) + 6 • [x 2 + y 1 + 16 x y) + 20 (y + x) (x 2 + y 2 + 80 x y) 
+ 70- (x 4 + y 4 +2b6x 3 y + 256xy 3 + 1296 x 2 y 2 + x 4 ) +••• (E.2) 

Note that all the coefficients of odd orders in x and y factor (x + y). 

The singular variety is the union of x y ■ (x — y) = together with the (x, y)- 
symmetric genus-zero algebraic curve which reads: 

2 s ■ (x - y) 4 - 2 8 ■ (x + y) ■ (x 2 + y 2 + 30xy) 

+ 2 5 - (3a; 2 + 3y 2 -62xy) - 2 4 • (x + y) +1 = 0. (E.3) 
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This genus-zero curve has the following polynomial parametrization: 

- ^ >= ^ 

In the y = x limit the singular variety (|E.3I) gives (1 — 64 a;) • (1 + 16 x) 2 = 0, in 
agreement with the singularities of the order-four Calabi-Yau operator (|E.1|) . 

Appendix E.2. Last Picard-Fuchs system with two variables 

A last example is the two- variables Picard-Fuchs system "above" the order-four Calabi- 
Yau operator (see ODE number 19 in [H])) 

529 9 4 - 23 x ■ (921 9 4 + 20 46 9 3 + 164 4 9 2 + 62 1 9 + 92) 

- x 2 • (380851 9 4 + 1328584 9 3 + 1772673 2 + 1033528 9 + 221168) 

- 2 x 3 ■ (475861 9 4 + 1310172 9 3 + 1028791 6 2 + 208932 9 - 27232) 

- 68 x 4 ■ (8873 9 4 + 14020 (9 3 + 5139 9 2 - 1664 6» - 976) 

+ 6936 x 5 ■ (3 9 + 4) ■ (3 9 + 2) ■ {9 + l) 2 . (E.5) 
The Picard-Fuchs system corresponds to the double series [35] 



EE 



(n + m)\ (2n + m)! (2m + n)\ 



n! 4 to! 4 



n— m— 

E(2to)! / r 1 m + 1, r \ 

■ 4 F 3 [[m + 1, m + -, 2m + 1, — ^— ], [1, 1, 1]; 4yj ■ x n 

m=0 

= 1 + 2 • (x + y) + (6 x 2 + 6 y 2 + 72 x y) + 20 • (x + y) ■ {x 2 + y 2 + 53 xy) 
+ 10- (1120a;?/ 3 + 1120a: 3 ?; +7x A + 7y 4 +4860a; 2 y 2 ) +••• (E.6) 

Note that all the coefficients of odd orders in x and y factor (x + y). 

The singular variety is the union of xy ■ (x + y) = together with the (a;, y)- 
symmetric genus-zero algebraic curve which reads: 

27- a; 2 ?/ 2 • {y + x) - [256 (a; 4 + y 4 ) +304a;?y • (a; 2 + y 2 ) + 69 a; 2 y 2 ] 
+ 8 ■ (y + x)- [32 (a; 2 + y 2 ) + 339 x y] 

- [96 (a: 2 +y 2 ) - 1261a:?/] +16- (y + x) -1 = 0. (E.7) 
with the simple rational parametrization (see section ([5])): 

ft 4 1 \ 

= \(t+l)(t + 2) (2t +1) 2 ' + (* + 2) 2 (2t (E,8) 

In the y — x limit, this singular variety reduces to (1— 54 x) ■ (1 + 11 a; — a: 2 ) 2 = 
in agreement with the singular points of (|E.5I) . 
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